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Abstract 

This paper explores the use of negative (i.e., repulsive) interactions in the abstract Tile As- 
sembly Model defined by Winfree. Winfree in his Ph.D. thesis postulated negative interactions 
to be physically plausible, and Reif, Sahu, and Yin studied them in the context of reversible 
attachment operations. We investigate the power of negative interactions with irreversible at- 
tachments, and we achieve two main results. Our first result is an impossibility theorem: after 
t steps of assembly, Sl(t) tiles will be forever bound to an assembly, unable to detach. Thus 
negative glue strengths do not afford unlimited power to reuse tiles. Our second result is a 
positive one: we construct a set of tiles that can simulate an s-space-bounded, i-time-bounded 
. Turing machine, while ensuring that no intermediate assembly grows larger than O(s), rather 

^"*^ | than 0(s ■ t) as required by the standard Turing machine simulation with tiles. 
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1 Introduction 

Tile-based self-assembly is a model of "algorithmic crystal growth" in which square "tiles" represent 
\& \ molecules that bind to each other via highly-specific bonds on their four sides, driven by random 

mixing in solution but constrained by the local binding rules of the tile bonds. Erik Winfree [11], 
based on experimental work of Seeman [8], modified Wang's mathematical model of tiling [10] 
to add a physically plausible mechanism for growth through time. Winfree defined a model of 
tile-based self-assembly known as the abstract Tile Assembly Model (aTAM). The fundamental 
components of this model are un-rotatable, but translatable square "tile types" whose sides are 
labeled with "glues" representing binding sites. Two tiles that are placed next to each other are 
attracted with strength determined by the glues where they abut, and in the aTAM, a tile binds 
^ to an assembly if it is attracted on all sides with total strength at least a certain threshold value 

r. 1 Assembly begins from a "seed" tile and progresses until no more tiles may attach. 

We study a variant of this model in which glue strengths are allowed to be negative as well as 
positive. This leads to the situation in which a stable assembly may become unstable through the 
addition of a tile that, while binding strongly enough to the assembly to remain attached itself, 
exerts a repulsive force on a neighboring tile, which is sufficiently strong to detach some portion 
of the assembly. This is formally modeled by allowing an assembly to break into two parts any 



*This research was supported in part by Natural Sciences and Engineering Research Council of Canada (NSERC) 
Discovery Grant R2824A01 and the Canada Research Chair Award in Biocomputing to Lila Kari. 

^University of Western Ontario, Dept. of Computer Science, London, Ontario, Canada, ddoty@csd.uwo.ca. 

* University of Western Ontario, Dept. of Computer Science, London, Ontario, Canada, lila@csd.uwo.ca. 

§ IRISA (INRIA), Campus de Beaulieu, 35042 Rennes Cedex, France, benoit.masson@irisa.fr. 

x The threshold r models the temperature at which insufficiently strong chemical bonds will break, such as those 
formed by Watson-Crick complementarity in DNA-based implementations of tiles. 



1 



time that the two parts have total connection strength less than r (i.e., if there is a cut of the 
interaction graph of strength less than r). Negative glue strengths were discussed as a plausible 
mechanism in Winfree's thesis [11], and explored theoretically in a more general model of graph- 
based self-assembly by Reif, Sahu and Yin [5]. We compare the results of [5] to the present paper 
in more detail later in this section. 

This paper has two main contributions, an impossibility result and a positive result. The 
impossibility result is that under the irreversible model, negative glue strengths are not sufficient 
to achieve perfect reuse of tiles as in [5]. It is tempting to believe that with negative glue strengths, 
the monotonic growth of the aTAM could be overcome to such a degree that a bounded set of tiles 
could be reused for arbitrarily long computations, 2 hence implementing the observation that "you 
can reuse space but you can't reuse time". Alas, you cannot reuse space (tiles) too much with 
irreversible reactions. In this section, we show that under the irreversible model of tile assembly, 
even with negative glue strengths, many tiles will be forever bound to an assembly, unable to detach. 
In fact, this number is linear in the number of assembly operations, so that after t operations, f2(i) 
tiles will be permanently bound to some assembly. 

The positive result is a construction attempting to make do with this limitation. For concrete- 
ness, our construction shows how to simulate a single-tape Turing machine. But the idea applies 
to the iterated computation of any function / that can be "computed with constant height" by 
a tile assembly system (a formal definition is given in Section 4). The function / mapping the 
configuration of a Turing machine to its next configuration is an example of one such function. 
Other examples include the incrementing or decrementing of a counter, or the selection of a uni- 
formly distributed random number from a finite set {1,2, ... ,n} using flips of a fair coin via von 
Neumann's rejection method, as shown in [4]. 

Our construction achieves the following property: if the Turing machine M being simulated on 
input x (with n = \x\) has space bound s(n) and time bound t(n), then 0(t(ri) ■ s(n)) tiles (meaning 
total count of tiles, which is greater than the number of unique tile types), mixed in solution, will 
simulate the computation of M on input x, and no intermediate assembly will grow to size larger 
than 0(s(n)). The impossibility result can be interpreted to imply that external energy must 
be supplied to break bonds between tiles if we wish to reuse them for computation. If we wish 
to limit the volume of a solution, and therefore the number of molecules it can contain (by the 
finite density constraint, see [9]) to 0(s(n)), then we cannot allow intermediate assemblies to grow 
larger than this value. Of course, by the impossibility result, many more than s(n) different such 
assemblies will form if t(n) » s(n) (for instance, when simulating a logarithmic-space, polynomial- 
time computation). With a mechanism to "vacuum" away junk assemblies and supply the external 
energy needed to break them up (a mechanism not modeled in the aTAM), these tiles could be 
reused, bringing down the required number of tiles from 0(t(n) ■ s(n)) to 0(s(n)). 

The main difference between [5] and the present paper is that [5] employs reversible reactions, 
and the present paper employs irreversible reactions. 3 Within the aTAM, the main difference 
between our model and [5] amounts to a difference in the definition of a legal attachment operation. 
In [5], the authors define a tile attachment to be legal if the tile attaches with strength r — 1 (in 
fact, they define it a bit differently but restricting attention to our construction and that of [5], 

2 subject, of course, to computational complexity constraints such as DTIME(i(n)) C DSPACE(2*'™'), based on the 
observation that configurations cannot repeat during the course of a halting computation 

3 [5] also uses a more general graph-based model of self-assembly, but this difference is less crucial than the 
reversibility issue. 
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this definition is equivalent). This is a phenomenon not modeled by the aTAM, but it is physically 
plausible to suppose that it occurs, though with less frequency than strength r attachments (see the 
kinetic TAM of [11]). Therefore the tile may detach after attaching since it is held with insufficient 
strength. But, if it first causes another tile or group of tiles to be bound with total strength less than 
r, then those tiles may also fall off, possibly resulting in stabilization of the original attachment. In 
the present paper, we define attachments to be legal only if they have strength at least r, whereas 
detachments may only happen between assemblies attached with strength at most r — 1. This 
difference implies that our impossibility result does not apply to [5], which can be considered an 
advantage of reversible interactions. But this advantage does not come without disadvantages: due 
to the second law of thermodynamics, their construction must necessarily be implemented as an 
unbiased random walk with equal rates of forward and reverse reaction, lest the entropy of the 
system increase with time if one direction is more favorable. Therefore their construction takes 
expected time n 2 to go forward n steps. 

We should also note that although [5] uses a more general model of graph-based self-assembly, 
this does not imply that their construction of an assembly system simulating a space-bounded 
Turing machine simulation is a stronger result than our construction. The more general model 
affords more power to aid in a construction, such as allowing non-planar interactions, in addition 
to the extra power of reversible interactions. Therefore, we emphasize that our positive construction 
is not merely a specialization of the construction of [5] to grid graphs. The construction of [5] is 
inherently non-planar and reversible. A major source of the effort in designing our construction 
was getting it to work in the plane and use irreversible attachments. Many similar (and simpler) 
constructions that superficially appear to do the same thing as our construction do not actually 
work in our model, as they introduce not only a desired cut of strength less than r, but also some 
undesired cuts of strength less than r, which, if detached, will ruin the construction. 

For color figures, see http://www.csd.uwo.ca/~ddoty/papers/niisa.pdf. 

2 Abstract Tile Assembly Model 

This section gives a brief definition of the abstract Tile Assembly Model (aTAM, [11]) with negative 
glue strengths. This not a tutorial on the aTAM; for readers unfamiliar with the model, please 
see [6] for an excellent introduction. Section A.l discusses the options we identified in defining an 
appropriate model and explains the choices we made. 

Z and Z + denote the set of integers and positive integers, respectively. Let G be a finite alphabet 
of glues. A tile type is a tuple t G G 4 , i.e., a unit square with a glue on each side. Associated with 
the tile types is a glue strength function str : G x G — > Z that indicates, given two glues g\ and 
g2, the strength str{g\,g2) with which they interact. We assume a finite set T of tile types, but 
an infinite number of copies of each tile type, each copy referred to as a tile. Let G(T) denote 
the set of all glues of tile types in T. An assembly (a.k.a., supertile) is a positioning of tiles on 
the integer lattice Z 2 (i.e., a partial function a : Z 2 T, where --■ » denotes that the function 
is partial). Each assembly induces a binding graph, a grid graph whose vertices are tiles, with an 
edge between two tiles if they are adjacent (i.e., are Euclidean distance 1 apart). 4 The assembly 
is T-stable, or simply stable if r is understood from context, if every cut of its binding graph has 

4 Previous papers model the binding graph as having edges only between tiles that interact with positive strength. 
In the present paper, the presence of negative glue strengths means that we must consider every possible interaction 
between adjacent tiles, whether positive, negative, or 0. 
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weight (strength) at least r, where the weight of an edge is the strength of the glue it represents. 
That is, the assembly is stable if at least energy r is required to separate the assembly into two 
parts. In this paper, where not stated otherwise, we assume that r = 2. 

A tile assembly system (TAS) is a 4-tuple T = (T,str,a,r), where T is a finite set of tile 
types, str : G(T) x G(T) — > Z is the glue strength function, a : 1? — + T is the finite and r-stable 
seed assembly, and r £ Z + is the temperature. Given a TAS T = (T, str,a,r), an assembly a is 
producible if either (base case) a = a, or (recursive case 1) a results from the r-stable attachment 
of a single tile to a producible assembly ("r-stable attachment" meaning that the cut separating 
the tile from the rest of the assembly has strength > r), or (recursive case 2) a consists of one side 
of a cut of strength < r of a producible assembly. Note in particular that a producible assembly 
need not be stable, but may be stabilized by attachments before it can break apart. An assembly a 
is terminal if a is r-stable and no tile can be r-stably attached to a. Let B C T be a set of "black" 
tile types. T is B-directed (a.k.a., B- deterministic, B-confluent) if it has exactly one terminal, 
producible assembly containing one or more tiles from B. 5 

To define reversible assembly at temperature r = 2 (as in [5]), it suffices to define attachment 
events with strength threshold r — 1 = 1, rather than strength threshold r = 2. This behavior is 
illustrated on Fig. 1(a), and can be compared with our new notion, whose evolution is shown on 
Fig. 1(b). 
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(a) Reversible model, as defined in [5]. 



(b) Irreversible model. 
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Figure 1: Two different implementations of negative interactions at temperature 2. The slanted bonds 
represent a strength of —1. In the reversible model, the tile t^ can attach with a total strength of 1 (one 
bond of strength 2 and one of strength —1) and hence is unstable, while with our definition, t' 3 is attached 
with a total strength of 2 and forces to to detach. 



3 Limitation of Tile Reuse with Irreversible Reactions 

If a is an assembly, define <&(«), the (negative) free energy of a, to be the sum of all glue strengths 
between adjacent tiles in the assembly. 6 In particular, an assembly consisting of a single tile has 
free energy 0. If S is a multiset of assemblies (such as that produced by a TAS with negative 
glue strengths, considering even the "junk" assemblies that are discarded after a cut), define the 
(negative) free energy of S to be the sum of the free energies of each assembly in S, denoted ${S). 

5 We define this notion of B-directedness but do not henceforth discuss it explicitly, since our construction simulates 
a general "computation", and B would depend on the goals of the computation being simulated. In our example 
construction in Section 4 of simulating steps of a Turing machine, B could, for instance, consist of the tile types that 
represent a halting state, so that only a terminal assembly representing the configuration of a halted Turing machine 
would be considered the result. 

6 The standard definition of free energy is the negative of this quantity, but as in [6] we use its negation so that 
the quantity will be positive for stable assemblies. Intuitively, it is the energy required to separate a into individual 
tiles, whereas the standard definition is the energy released by such a separation. 
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Note that even postulating an infinite count of tiles, after a finite number of operations, only finitely 
many assemblies in S consist of more than one tile, and each of these is a finite assembly. Therefore 
$>(S) < oo for any multiset S of assemblies producible by a TAS, even in the case that \S\ = oo 
(such as the initial multiset consisting of a countably infinite number of copies of each individual 
tile type). 

When we discuss the "number of steps" for the assembly process of a TAS, we mean the total 
number of attachment and detachment operations that have been applied so far. We do not claim 
that this is a proper model of "running time", but it is convenient to think of attachment and 
detachment events as discrete and equally-spaced steps, even though they may happen in parallel 
or with interval times governed by a continuous distribution. 

Theorem 3.1. Let T be a TAS, and let S be a multiset of assemblies producible by T after t £ N 
steps. Then $(S) > t/2. 

A proof of Theorem 3.1 is given in Section A. 3. 

In particular, since the glue strengths are bounded above by some constant s, this implies that 
after t steps, at least t/(2s) sides of tiles are bound. With the finite tile count assumption, once t 
is sufficiently large that t/(2s) exceeds the total number of sides available (i.e., 4 times the total 
number of tiles in solution), no more sides are available for binding, and self-assembly grinds to a 
halt. This is the sense in which a finite number of tiles cannot be reused indefinitely. 

Some seemingly straightforward attempts to prove Theorem 3.1 fail in ways that illustrate 
potentially nonintuitive properties of negative glue strengths. It is not true, for instance, that the 
free energy increases monotonically, since it drops whenever a cut of positive strength is detached, 
so a straightforward inductive argument fails. Furthermore, it is not even true that the free energy 
decreases by at most a constant in between consecutive periods of increase. Even with fixed glue 
strengths (4 suffices), for each n G N, it is possible to construct a tile set with the property that there 
are two states of the system Si and Sj, with state Si preceding Sj, such that $(Si) > $(Sj) + n. 
But, Theorem 3.1 implies that any attempt to create such a cascade of detachments that drops 
the free energy by n requires first attaching even stronger - and ultimately unbreakable - bonds 
required to set up the state Si. That is, the free energy can fall arbitrarily far, but in order to do so 
it must first climb more than twice as high as it will eventually fall. This phenomenon is illustrated 
in Fig. 2, where the last stable addition of a tile leads to an arbitrary decrease of the free energy. 
This was made possible by the use of stronger strength-4 bonds prior to this event. 

There is a natural thermodynamic interpretation of Theorem 3.1: work done by tiles on tiles, in 
an irreversible manner, increases the entropy of the system by the second law of thermodynamics, 
thus decreasing the potential energy available to do more work. Therefore, any potential energy 
stored in the unattached glues is eventually permanently used up if external energy is not supplied 
to break these bonds. In our main construction, many junk assemblies are created that are no 
longer useful once the tiles in them have been used once. Theorem 3.1 tells us that no amount of 
cleverness will allow us to break up those assemblies and reuse the tiles solely through design of 
tile types with negative glues; some external force must be supplied to break them apart using a 
mechanism not modeled in the aTAM. 

Of course, Theorem 3.1, interpreted in light of the molecular interactions that are being modeled 
by the aTAM, should not be surprising to any physicist. But we believe it is important to formally 
establish the truth of such a statement within the model. One develops more confidence in a 
model of reality when it tells us something already known about reality (e.g., the Positive Mass 
Theorem [7]). 
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n tiles n tiles n tiles 




(a) Starting from the seed a, n — 1 (b) From the gray tile, n — 1 light (c) The initial tiles are removed by neg- 

tiles attach on the left, before gray tiles attach on the right, ative bonds of strength —2, and re- 

a gray tile can be attached at using the bonds on the top placed by n new tiles until t is at- 

the bottom. and on the left. The assem- tached. The next step is the detach- 

bly also goes around these new ment of the gray tile, followed by the 

tiles counter-clockwise, until a detachment of the n — 1 light gray 

is reached. tiles one after the other. 



Figure 2: A possible evolution where the stable addition of one tile (marked t) can lead to n tiles detaching 
one after the other (black arrows), hence reducing the free energy by n. As usual, the strength of bonds is 
represented by the number of segments between tiles, slanted bonds indicating a negative strength. Note 
that this is only one among many possible evolutions, since there may be several cuts of strength lower than 
t that can be removed. In particular, at the last step, all the gray and light gray tiles can detach as one 
unique big cut, which will in turn break into pieces. 



Theorem 3.1 does not apply to the negative glue strength construction of Reif, Sahu, and Yin [5], 
because their model allows reversible reactions. Intuitively, the reversibility of reactions implies 
that attachment and detachment have symmetric effects on the free energy. But this also implies 
that their system requires driving the system forward through an unbiased random walk, taking n 2 
steps on average to proceed by n net forward steps. Any attempt to speed up the reaction to make 
the forward rate of reaction faster than the reverse rate of reaction would introduce the imbalance 
in their respective effects on free energy that allows our proof to work. Therefore this tradeoff in 
speed versus reusability of tiles is fundamental. 



4 Turing Machine Simulation 

Throughout this section, fix some finite alphabet S. We first describe the class of functions that 
we will compute, which are intuitively those computable by a constant number of rows of assembly 
(although the number of columns is unbounded) in the standard aTAM. See [4], for example, for 
a formal definition of the standard aTAM model. Briefly, it is the same as the model defined in 
Section 2, but glue strengths are non-negative and are only positive between equal glues. 

Informally (see Section A. 4 for a formal definition), a function / : S* — > S* is constant-row 
computable if there is a tile set T and a constant c € N such that, for each i G S*, if a horizontal 
row a x of tiles from T that "represents" a: in a straightforward way is used as the seed assembly, 
then the resulting tile assembly system, at temperature 2, will grow into an assembly whose top 
row "represents" f(x). Furthermore, the entire assembly will have height exactly c, and the left 
most tile of each row will have the same horizontal position as the left-most tile of the bottom row. 

The widths of the rows representing the input and output may be different (i.e., possibly 
\x\ |/(x)|). In this case, we require only that the leftmost and rightmost tiles of each row have 
their glues specially marked to distinguish them from the tiles interior to the row. 
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Figure 3: High-level overview of assembly for computation of a constant-row computable function /. 



Our construction shows how to design a tile set that will compute iterations of any constant- 
row computable function /, ensuring that no intermediate assembly grows larger than the size 
of the input or output processed by any individual invocation of /. Examples of such functions 
include the function / that, given a configuration of a single-tape Turing machine outputs the next 
configuration of this Turing machine, or that increments a counter represented in binary. 

Figure 3 shows a high-level overview of the entire construction, in terms of a general constant- 
row computable function /. For concreteness, think of / as the function that, given a configuration 
of a t-time-bounded, s-space-bounded, single-tape Turing machine, outputs the next configuration 
of this Turing machine (extending the tape on the right side only). The construction proceeds as 
follows, each label corresponds to a picture in Figure 3. 

(a) First, the scaffold tiles (green) connect to the x data assembly (white). The scaffold tiles 
initiate the computation of / (blue). 

(b) The scaffold "detects" when the computation is finished, in the sense that the green row 
above f(x) tiles cannot complete until all of f{x) is present. Then the scaffold tiles grow back 
to the first scaffold tile to initiate the removal of f(x) from the tiles surrounding f(x). 

(c) The removal tiles (red) each use a negative glue strength against the tile "in front of" (on 
the path show by the arrows) it, and once this tile is removed, a new removal tile grows in 
its place to continue the removal. The path and bond placements and strengths are carefully 
chosen to ensure that no portion of f(x) is removed, until the last step when f(x) detaches 
whole from the rest of the tiles. 

Note that since / is constant-row computable, the height of the scaffold and removal parts are 
bounded by a constant and therefore may be hard-coded into the tile set, whereas special glues 
mark the horizontal endpoints so that the length of x and f{x) are not constrained. 

The simulation of the Turing machine for t steps will then consist of executing this assembly 
process for t iterations, using the output assembly f{x) as the input assembly x for the next 
execution. After each iteration, the width of the remaining "junk" assembly is a constant plus 
O(l) + max{|x|, |/(x)|}, and the height is constant since / is constant-row computable, so the size 
of the intermediate assemblies is 0(max{|rc|, |/(x)|}). 
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Figure 4: Example of tiles implementing the computation step. Arrows within tiles show order of growth. 
In this case / is constant-row computable with constant c = 1. The first and last copy rows, shown in lighter 
shade than the center computation tiles, are always present no matter the function /, and their placement is 
initiated by the scaffold tiles. However, there is no interaction between the center computation and scaffold 
tiles. Note that the data tiles are two rows with strength I glues; this is to make them stable at temperature 
2 but not producible (without additional scaffolding) as they would be if they were a single row connected 
with strength 2 glues. 



Figures 4, 5, and 6 give some more details for the three main steps of Figure 3, respectively 
(a), (b), and (c), using the specific example of / mapping a configuration of a single-tape Turing 
machine to its next configuration. 

Figure 4 shows an example of tiles implementing step (a) of Figure 3, i.e., the computation 
of /. The example shows one transition of a single-tape Turing machine, with tape contents 01_0 
(_ standing for blank), in state q, with tape head on the rightmost cell, transitioning to state p, 
moving the tape head right, changing the cell's symbol from to 1, and encountering a blank on the 
new rightmost cell. In this case, a new rightmost cell is needed, illustrating how our construction 
handles dynamically changing space requirements, but if the tape head were further left in the row, 
it would simply fill in copy tiles to the right, just as to the left as shown above, and the row would 
stay the same width. At the start and end of a computation, the configuration is copied so that any 
strength > 1 bonds used in the computation are on the interior of the computation tiles, ensuring 
that only strength-1 bonds must later be broken to separate the data tiles. Each data assembly on 
either end of the computation tiles is represented by a two-row assembly with only single-strength 
bonds on its interior, which ensures that when detached, the data assembly will be stable, but that 
it cannot form on its own without help from the scaffold tiles (which would happen if it were only 
a single row connected with strength-2 bonds). Each vertical position is hard-coded into the tile 
set; i.e., the scaffold tile set "knows" the required height to compute /. However, the absolute 
horizontal positions are not encoded into the tiles, only the leftmost and rightmost tiles of the 
configuration are specially marked, and all interior tiles representing the same data are identical. 

Figure 5 shows the tiles implementing step (b) of Figure 3, positioning the tiles for cleanup. The 
top two rows must use cooperation to tell where the end of the row underneath is, since the width 
of the output row is unknown. The strengths of bonds on the leftmost downward-growing column 
must be sufficiently large to ensure that only the proper cut is made when the first negative-strength 
glue is applied. 
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Figure 5: Tiles that position the cleanup tiles. Here the "copy" tiles from Figure 4 are depicted in the same 
shade as the computation tiles; now that fix) has been computed our goal is to remove all of them from 
the subassembly representing fix). The order of growth of the scaffold tiles ensures that cleanup does not 
begin until all of fix) is present. 
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Figure 6: Tiles that "clean up" the connections between the output data and the scaffold and computation 
tiles to separate them and allow the data tiles to be computed on again. 
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Figure 6 shows the tiles implementing step (c) of Figure 3, "cleaning up" by removing the output 
f{x) from the scaffold, computation, and x tiles. Though not shown, negative strength interactions 
are necessary between the second-to-top row of computation tiles and some of the right-growing 
cleanup tiles, to ensure that the right end of the row is properly detected. That is, there are two 
types of cleanup tiles growing right, one to detach the interior tiles, and one to detach the final 
rightmost computation tile. Since the east-west bonds between cleanup tiles are greater than 1, the 
negative north-south glue strengths between interior cleanup tiles and the second-to-rightmost blue 
tile - and between the rightmost cleanup tile and the interior computation tiles - must be strength 
-2 to ensure that the second-to-rightmost blue tile cannot stably attach except where intended. 

5 Conclusion 

There are other uses of negative glues in the aTAM. For instance, we are able to improve the best 
known tile complexity (number of tile types) required to uniquely assemble a "thin" rectangle, 
i.e., fflinxfc rectangle with k < log nj (log log n — log log log n). In the standard aTAM the tile 

l/fc 

complexity of this shape is known to be f2(^V— ) [1]- With the model of negative glue strengths we 
are able to improve this to O (\/logn) , by first building a thick rectangle and using negative glues 
to "cut out" a thinner rectangle of the same length. We refer to Section A. 5 for more applications. 

Other questions related to this work include the experimental aspects of such a model, for 
example, how repulsive forces can be realized on DNA tiles, and how to "clean" and "recycle" the 
junk introduced during the assembly. 

Acknowledgments. The authors are grateful to Shinnosuke Seki for insightful discussions. 
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A Appendix 



A.l Issues with Choice of Model 

There are many variations on the model of tile self-assembly with negative glue strengths. We 
identify six (somewhat) independent binary choices to be made: 

seeded/unseeded: In the seeded model, assembly starts from a specially designated seed tile. In 
the unseeded model, assembly may start with any tile. 

single-tile addition/two- handed assembly: In the single-tile addition model, tiles are added 
one at a time to an assembly. In two-handed assembly [1], two assemblies, each possibly 
consisting of multiple tiles, may attach to each other. 

irreversible/reversible: In the irreversible model, for a tile/assembly to attach to an assembly, 
it must bind with strength > r (though it may cause another cut of the resulting assembly to 
have strength < r). In the reversible model (used by Reif, Sahu, and Yin [5]), a tile/assembly 
may attach with strength < r, implying that it may detach (reverse), but may also cause 
another cut to detach. In [5], most attachment events have the property that they cause 
precisely two cuts to have strength r — 1, that of the attachment event itself, and another 
desired cut, and if the latter cut detaches, then the former cut now has strength r. Hence, 
on the assumption that each cut is equally likely, assembly proceeds in a unbiased random 
walk (thus taking expected time n 2 to go forward n steps). 

detachment precedes attachment/detachment and attachment in arbitrary order: If de- 
tachment precedes attachment, this means we assume that whenever a negative-strength glue 
creates an unstable cut, every detachment event that can occur, will occur, until we are left 
with nothing but stable assemblies, before the next attachment event occurs. The more re- 
alistic model, detachment and attachment in arbitrary order, assumes that the operations of 
attach and detach are both legal at any stage. Therefore, if a negative-strength glue creates 
an unstable cut, but this cut could be stabilized if an attachment secures it in place before it 
can detach, then we assume that this could happen. 

finite tile counts/infinite tile counts: Most models of tile self-assembly assume an infinite num- 
ber of tile types, since a very large number can be easily created in a short time. However, 
if we wish to use negative glue strengths to implement space-bounded computation, and 
the finite density constraint (see [9]) implies that the solution volume must be proportional 
to the number of molecules it contains, then we must constrain the number of molecules. 
See the next point for a discussion of another practical issue with using bounded space for 
computation with tile assembly. 

tagged result /tagged junk: This choice relates to how we designate what is the "result" of 
assembly. In the tagged result model, we designate a subset of tile types to be "black" , and 
state that a result assembly is any terminal assembly with a black tile in it. This allows us 
to separate the junk from the result after assembly is complete, but does not allow junk to 
be removed during assembly, since the black tile may not be attached to the result assembly 
until the very end of the assembly process. In the tagged junk model, we say that any 
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assembly (terminal or not) with a black tile is junk and may be removed. 7 Furthermore, 
every producible assembly must have the property that it may either grow into the result, 
or will grow into an assembly tagged with a black tile as junk. That is, at any point in the 
assembly process, we could "vacuum" away all assemblies tagged with a black tile, knowing 
that anything remaining is not (yet) junk and can be re-used. This allows a computation 
with space bound s and running time t > s to be done in volume O(s), so long as one uses 
detachment to ensure that no intermediate assembly grows larger than O(s), by supplying 
just enough copies of tiles that are not swept away, and periodically sweeping away the tagged 
junk while supplying fresh unattached tiles to carry out more work. 

Note that neither of these tagging models are typically used in other tile self-assembly papers, 
but the special implications of negative glue strengths imply that we cannot simply follow 
the convention that the seed identifies the result, as discussed below. 

Of these, the first three choices are incomparable in terms of power: each choice affords both 
advantages and disadvantages in terms of designing a tile assembly system. The latter three choices 
are more clear: for each choice, one option makes implementing a correct design strictly more 
difficult, but results in a more robust and realistically implementable construction. Respectively, 
these choices are 1) detachment and attachment in arbitrary order, 2) finite tile counts, and 3) 
tagged junk 

These choices are not completely independent. For instance, there are two (mathematical) uses 
for a seed tile: 1) to identify legal attachment events in single-tile addition: attachment is legal if it 
is between a single tile and an assembly containing the seed, and 2) to identify the result of assembly 
with two-handed growth: one must allow attachment of tiles separate from the seed, but identifies 
legal results as terminal assemblies containing the seed. We separate out these uses by making the 
method of "tagging" the result independent of the seed tile, using the seed only to define what 
assemblies are defined as producible. With negative glues implying the ability to remove the seed 
from an assembly, we must take care to allow attachment events in the seeded model not only with 
assemblies containing the seed, but with those derived from the seed. 

For this extended abstract, we use the model of single-tile addition, irreversible, seeded, de- 
tachment and attachment in arbitrary order, infinite tile counts, and tagged result. We believe it 
is possible to modify our construction to allow unseeded, two-handed assembly, finite tile counts, 
and tagged junk (the latter two allowing a more realistic implementation of "space-bounded com- 
putation" as discussed above), but we delay this construction until the full version of this paper. 

A. 2 Definitions of Other Models 

To define unseeded assembly, it suffices to drop a from the definition of TAS, and define the base 
case of a producible assembly as any individual tile. To define two-handed assembly (a.k.a., multiple 
tile model [1]), it suffices to change the first recursive case to state that legal attachment events are 
between any two producible assemblies, such that they can be positioned in such a way that the cut 
separating them has strength > r (i.e., can be stably attached). Then, an assembly a is terminal 

7 Actually it is more realistic to assume something like two different black tiles, placed at a certain recognizable 
orientation (such as next to each other) to enforce that individual black tiles in solution are not interpreted as 
identifying junk until they actually attach to something, but we do not dwell on this issue since we use the tagged 
result model. 
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if for every producible assembly /?, a and (3 cannot be stably attached. Figure 7 illustrates the new 
behaviors allowed by the two-handed variant. 



to 





■*3 


to- 


-H- 



(a) 



(b) 



(c) 



(d) 



Figure 7: Typical example of two-handed assembly, at temperature r = 2. The segments between tiles 
represent the bonds, the number of segments encodes the strength of the bond (here, 1 or 2). In the seeded, 
single tile model with seed a = to, the assembly at step (b) would be terminal. 



A. 3 Proof of Theorem 3.1 

Theorem 3.1. Let T be a TAS, and let S be a multiset of assemblies producible by T after t £ N 
steps. Then $(5) > t/2. 

Proof. Suppose that T has a seed of size 1; otherwise, the free energy we derive for step t will 
be even higher, so this assumption does not harm the proof. For i G {0, 1, . . . , t}, let Si denote 
the multiset of assemblies after the first i operations, so that St = S and So is the multiset of 
individual unattached tiles. Note that $(So) = 0. Let A C {1, 2, ... ,t} be the indices of attachment 
operations, and let D = {1, 2, . . . , t} — A be the indices of detachment operations, so that operation 
i changes Si-i to Si. 

Each attachment operation increases the free energy by at least r for a system operating at 
temperature r, since we require a tile/supertile attachment to have the property that the cut 
between it and the rest of the assembly has strength at least r, and the edges of this cut previously 
each contributed to the free energy since they were all unbound. So for Sj_i leading to Si via 
attachment, — $>(Si-i) > r. For each detachment operation, the greatest strength cut that 

could be broken to create the detachment has strength r — 1; stronger cuts cannot be broken. This 
implies the free energy decreases by at most r — 1 during a single detachment operation. 8 So for 
Si-i leading to Si via detachment, — 3>(Si_i) > — (r — 1). Amortizing over all operations, 9 

we have that 



> £^ + £^-(^-1) = |A|r-|Z>|(r-l). 

Although we posit an infinite number of copies of each tile type, during the first t steps at most 
t+1 tiles, denoted as the multiset S' , can actually participate in assembly operations. Let q denote 
the total number of assemblies in Si that consist of tiles from S' , so that cq is simply the total 
number of initial tiles that will participate in the first t steps. Each attachment event decreases Cj 
by one, and each detachment event increases q by one. Since q < c$ for all i (we cannot have more 

8 A single attachment of a tile with negative glue strength can potentially cause a cascade of detachments that, 
put together, lead to a large decrease in free energy. However, these are each considered separate detachment events. 

9 See [2, Section 17.3] for a discussion of amortized analysis, which is a fancy phrase for writing the following sum 
in this form. 
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assemblies than there are tiles), this implies that for all i, the number of attachment events in the 
first i steps is at least the number of detachment events in the first i steps. Therefore \A\ > \D\ 
and since |^4| + \D\ = t, we conclude \A\ > t/2, whence the above inequality tells us that 



The proof works for two-handed systems, and since single-tile addition systems are simply 
two-handed systems with an extra constraint on legal attachment operations, the proof applies to 
single-tile addition systems as well. The proof also works both for seeded and unseeded systems. 
As there is no property of the 2D plane or grid graphs employed in our proof, the proof applies to 
the irreversible version of the graph-based self-assembly model studied in a reversible context by 
Reif, Sahu, and Yin [5] . 

A. 4 Formal Definition of Constant-Row Computability 

The following is a formal definition of the class of functions that we call constant-row computable. 
The definition employs the standard model of nonnegative glue strengths, in which strengths are 
only positive between equal glues, so the glues themselves have both a "type" and a strength. 

Definition A.l. Let T be a set of tile types, and let e : T — > S. We say that a row of tiles 
(a connected subassembly of some assembly of height 1) ti, <2, ■ ■ ■ , tfc e-encodes a string x € S fc 
if e(ti) = x[l],e(t2) = x[2], . . . ,e(tk) = x[k], where x[i] G S is the i th symbol in x. A function 
/ : S* — > S* is constant-row computable if there exists a tile set T, a function e : T — > S, and a 
constant c such that, for each x £ £*, there is a height- 1 stable assembly o~ x : 1? — > T e-encoding 
x such that the tile assembly system T = (T, str,o~ x ,2) (with str(g\,g2) > -4=>- g\ = 52) has 
the unique terminal assembly a, the height of a is c, the bottom row of a is a x , the top row of a 
e-encodes f(x), and the leftmost column of any row of a is no further left than the bottom row. 

A. 5 Other Uses of Negative Glue Strengths 

Space-bounded "computation" as an end goal is not the only application of negative glue strengths, 
of course. Doty, Lutz, Patitz, Summers, and Woods [4] study the problem of generating uniform 
random distributions on the finite sets using the independent flips of a fair coin afforded by the 
random selection of competing tile types in the aTAM (a non-trivial problem when the cardinality 
of the set is not a power of the number of competing tile types), and find a tradeoff between the 
closeness to uniformity of the distribution obtained and the space required for sampling. They 
exhibit a construction imposing a perfectly uniformly distribution on the set {0, 1, . . . , n — 1} that 
assembles a structure of width [lognj + 1 and expected height at most 2, essentially implementing 
von Neumann's rejection method of flipping [log n\ + 1 fair coins repeatedly and stopping the first 
time that they encode a number smaller than n. It is very unlikely (probability at most 2~ 20 ) 
to take more than 20 attempts. But using this method in a construction such as that of [3], 10 in 

10 The main construction of [3] shows how a "universal" tile set can be constructed that can be "programmed" 
through appropriate selection of a seed assembly to simulate the growth of any tile assembly system in a wide class of 
systems termed "locally consistent" (see [3] for details). In this discussion, we are concerned only with the fact that 
the construction of [3] 1) requires random numbers to be generated in a bounded space at many points throughout 
assembly, and 2) would be improved if the distribution of these numbers were perfectly uniform instead of "close to 
uniform" as in [3]. 



$(S) > \A\t - \D\(t - 1) > \A\t - \A\(t - 1) = \A\ > t/2. 



□ 
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which many (perhaps more than 2 20 ) copies of this experiment repeat throughout assembly, could 
increase the likelihood of growing too high. Even a single occurrence of a too-high subassembly 
will destroy the entire construction. Though we omit the details in the extended abstract, it is 
straightforward to augment the construction of [3] (which uses a variant of the random number 
selector of [4]) with negative glue strengths to implement perfectly uniform selection of random 
numbers, thus improving the fidelity of the simulation of [3], while providing an absolute guarantee 
on the space bound. 
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